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摘 要: 基于 无 人 机 可 见 光影 像 对 古 尔 班 通 古 特 沙漠 地 表 类 型 信息 进行 提取 ,运用 面向 对 象 的 多 尺度 分 割 , 在 提取 
样本 的 光谱 形状、 纹理 、 植 被 指数 特征 的 基础 上 ,建立 规则 提取 地 表 类 型 信息 。 结 果 表 明 : 荒漠 地 表 类 型 不 同 ， 
最 佳 分 割 尺 度 不 同 ;@ 不 同 荒漠 化 程度 地 表 类 型 特征 相似 ,无 法 运用 单个 特征 进行 区 分 , 需 选 用 多 种 特征 组 合 提取 
地 表 类 型 ;(3) 面向 对 象 的 多 尺度 分 割 方法 相对 于 基于 像 元 的 最 大 似 然 法 分 类 有 明显 提高 ,面向 对 象 轻 度 沙漠 区 总 
体 分 类 精度 为 93.00% ,中 度 沙漠 化 区 为 91.83% ,重度 沙漠 化 区 为 93.50% , 较 基于 像 元 的 最 大 似 然 方法 分 别提 高 


了 10.3496 ,11. 8696 和 12.50% 。 表 明 针 对 无 人 机 可 见 光影 像 ， 


押 向 对 象 的 多 尺度 分 割 方法 能 高 精度 地 提取 蕊 漠 地 


表 类 型 信息 。 
关键 词 : 无 人 机 影像 ; icis: 地 


表 类 型 ; 

荒漠 化 是 一 个 全 球 性 的 环境 问题 ,对 荒漠 化 遥 
感 监测 始终 是 国际 上 关注 的 一 个 热点 y EE 
漠 化 主要 发 生 在 内 蒙古 和 西北 地 区 , 虽 有 一 定 的 组 
解 ,但 仍然 没有 得 到 根本 性 的 遏制 ”) 。 为 了 能 够 快 
速 获取 荒漠 化 趋势 ,通常 应 用 遥感 技术 , 因 它 监测 范 
围 大 ,时 间 周 期 短 ,能 够 快速 获取 地 面 信 息 的 特点 ， 
可 以 提高 对 荒漠 化 监测 的 效率 所。 但 随 着 荒漠 化 
研究 的 深入 ,发 现 作 为 荒漠 化 判 据 的 植被 盖 度 程度 
毕 合 指标 不 仅 受到 关注 ,而且 荒漠 生态 系统 的 生物 
多 样 性 地衣 等 空间 分 布 特征 也 日 益 得 到 荒漠 人 研究 
者 的 关注 5 。 目 前 ,已 经 有 很 多 传统 传感器 数据 
(如 MODIS Landsat 数据 等 ) 应 用 于 车 漠 化 监测 ,但 
它们 仅 能 捕 提 像 元 内 综合 的 植被 覆盖 信息 ,而 对 
于 荒漠 化 更 为 具体 的 信息 却 无 能 为 力 。 无 人 机 的 发 
展 改 变 了 这 一 难题 , 相 比 较 卫 星 遥 感 ,无 人 机 影像 有 
着 丰富 的 几何 纹理 信息 ,具有 捕捉 细节 信息 的 
能 力 。 

传统 的 分 类 方法 主要 是 基于 像 元 的 遥感 解 译 技 
术 , 多 应 用 于 中 低 分 辨 率 的 影像 ,并 取得 较 好 的 效 
果 , 而 在 高 分 辨 影像 中 容易 产生 椒盐 效应 ”") 。 针 对 
高 分 辩 率 影像 特点 提出 一 种 面向 对 象 分 类 的 方法 ， 
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它 能 捕 提 地 表 的 纹理 .几何 .颜色 等 信息 局 。Ke 
等 中 运用 面向 对 象 的 方法 ,对 黑 堡 纪念 森林 和 吡 邻 
的 纽约 中 部 林地 进行 分 类 ,结果 表明 面向 对 象 方法 
分 类 精度 高 于 基于 像 元 的 分 类 精度 。Myint 45079 基 
于 像 元 和 面向 对 象 方法 ,对 亚利桑那 州 风 凰 城 的 中 
心地 区 进行 分 类 ,结果 表明 面向 对 象 比 基于 像 元 的 
分 类 方法 好 。 这 些 研究 表明 ,面向 对 象 分 类 方法 要 
优 于 基于 像 元 分 类 方法 。 目 前 ,无 人 机 作为 新 兴 
术 在 很 多 领域 都 得 到 了 应 用 ,如 水 污染 监测 、 草 
原 分 类 调查 :2 、 森 林 资 源 "” 与 农作物 识别 "等 ， 
在 细节 信息 获取 上 展现 出 一 定 优 势 。 在 土地 利用 方 
面 ,以 无 人 机 影像 为 数据 源 ,运用 面向 对 象 的 方法 促 
使 分 类 精度 提高 到 91.30%"]。 在 林地 方面 ,利用 
无 人 机 的 优势 ,用 面向 对 象 的 方法 提取 植被 信息 , 精 
度 达到 90. 00% 以 上 , 优势 突出 5 。 在 农田 方面 ， 
无 人 机 影像 作为 数据 源 ,结合 面向 对 象 的 方法 ,能 够 
高 效 精 确 地 识别 农田 '""。 但 在 荒漠 化 方面 ,国内 外 
却 鲜 见报 道 ,尤其 是 在 国内 。 

古 尔 班 通 古 特 沙漠 属 典型 的 温带 荒漠 ,降水 稀 
少 ,气候 干 燥 , 革 发 强烈 避 ] 。 在 特殊 气候 影响 下 ,地 
表 植 被 盖 度 达到 30. 00% 左右 ,植被 不 仅 有 灌木 梭 
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梭 和 白 梭 梭 等 "”, 林 下 形成 了 片 状 的 草本 ,部 分 裸 
Vo e DE JY JO C Mo UU ,荒漠 地 表 形 成 高 空间 
异 质 性 特征 。 如 何 获取 鞠 漠 地 表 空 间 异 质 性 特征 是 
研究 车 漠 化 空间 演变 的 关键 环节 。 本 文 以 古 尔 班 通 
古 特 沙漠 为 研究 对 象 , 以 无 人 机 影像 为 数据 源 ,运用 
面向 对 象 的 分 类 方法 对 古 尔 班 通 古 特 沙漠 地 表 类 型 
进行 提取 ,为 荒漠 化 空间 演变 研究 提供 方法 参考 。 


1 材料 与 方法 


1.1 研究 区 概况 

古 尔 班 通 古 特 沙漠 地 处 准噶尔 盆地 , 玛 纳 斯 河 
以 东 及 乌 伦 古河 以 南 , 天 山 经 济 带 北 葛 , 位 于 
44?11' ~46°20'N, 84°31’ ~ 90?00' E 之 间 , 面积 约 
4.88 x 10* km ,是 我 国 面积 最 大 的 固定 和 半 固 定 沙 
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图 1 研究 


Fig. 1 Geographical location of the study area 


1.2 数据 来 源 和 处 理 

无 人 机 可 见 光 影像 获取 的 时 间 为 2013 年 7 H 
18 日 ,天气 晴朗 ,无 风 无 云 ,选择 正午 ,以 避免 太阳 高 
度 角 所 引起 的 阴影 对 影像 质量 造成 的 影响 。 航 高 为 
500 m , 预 设 航向 重 三 度 为 65% , 257 [n] EAEN 40% , 
相机 曝光 时 间 为 1/1 600 s ,型 号 为 SONY NEX -7, 在 
设 定 航 高 下 影像 空间 分 辩 率 为 0. 15 m。 无 人 机 获 
取 可 见 光影 像 , 其 中 包括 红 、 绿 、 蓝 3 个 波段 。 对 航 
摄影 像 进行 拼接 .几何 校正 及 影像 裁剪 等 ,结合 航拍 
无 人 机 GPS 记录 坐标 等 参数 ,用 SEM ( structure from 


ChinaXiv 合 作 期 刊 


漠 。 年 降水 量 70 ~ 150 mm ,蒸发 量 在 2 000 mm 以 
上 ,年 平均 气温 6 ~ 10 % 。 土 壤 以 固定 FEE KY 
土 占 绝对 优势 。 灌 木 植被 主要 有 梭 梭 (Haloxylon 
ammodendron ) 、 蛇 麻黄 ( Ephedra distachys ) | b 19 2& 


( Calligonum mongolicum ) , ÉE Œ #4% ( Reaumuria. son- 


garica) 等 。 短 命 和 类 短命 植物 发 育 良好 , 主要 植物 
PARKA F JL (Erodium oxyrrhynchum ) PER 
草 ( Carex physodes) 、 粗 柄 独 尾 草 ( Eremurus inderien- 
sis) fA RŽ ( Ceratocarpus arenarius) C? 等。 在 沙漠 
南 缘 到 腹地 设置 3 个 样 区 ,分 别 位 于 研究 区 南部 .中 
部 和 北部 ,实地 选择 典型 样 地 :GD 轻 度 沙漠 化 (固定 
沙丘 ) ;@@ 中 度 沙漠 化 ( 半 固 定 沙 丘 ) ;@@ 重度 沙漠 
ERD E) 。 

遥感 解 译 样 区 :3 个 样 区 ,每 个 样 区 1. 5 km x 
1.5 km。 如 图 1 所 示 : 


区 地 理 位 置 


"s 


motion ) HE 对 影像 进行 拼接 ,采用 多 项 式 几 何 
校正 。 选 用 最 邻近 像 元 法 进行 影像 重 采样 ,平均 误 
差 控 制 在 0.5 个 像 元 以 内 ,精度 可 靠 , 满足 研究 
2013 年 7 月 18—20 日 获取 野外 采样 点 数据 ， 
根据 研究 区 地 物 分 布 ,每 个 样 区 布设 了 10 个 涵盖 研 
究 区 典型 地 物 的 20 m x20 m 样 方 , 样 方 记录 内 容 包 
括 中 心 点 GPS 坐标 .地 物 类 型 灌木 冠 幅 等 。 
1.3 研究 方法 
采用 面向 对 象 的 多 尺度 分 割 方法 ,在 分 类 体系 
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的 基础 上 ,对 尺度 分 割 参数 特征 进行 分 析 ; 基 于 样本 
的 光谱 形状、 纹理 、 植 被 指数 特征 ,分 析 荒 漠 地 表 的 


查 ,将 研究 区 地 表 类 型 分 为 4 类 : 裸 阔 灌木. 结 皮 、 
Ric (3 1), 


可 见 光影 像 特征 ,并 建立 规则 ;采用 面向 对 象 的 多 尺 
度 分 割 方法 与 基于 像 元 的 最 大 似 然 法 分 别提 取 芳 浇 
地 表 类 型 信息 ,应 用 精度 评价 方法 对 两 种 分 类 方法 
结果 进行 对 比分 析 。 研 究 路 线 如 图 2 所 示 。 


上 面向 对 象 信息 提取 方法 
" _ | 最 优 分 割 | | 影像 | | 建立 规则 | | 分 类 | | 
T T Reta an | “| 提取 信息 | “| 结果 | ;| | 结果 
机 Ly 和 | 2 
M || 基于 像 元 信息 提取 方法 评价 


| Lm 分 类 | ; 
Pra muel JRA 


图 2 研究 路 线 图 


Fig.2 Diagram of the research route 


1.3.1 多 尺度 影像 分 割 ” 影 像 分 割 采用 多 尺度 分 
割 ,使 同一 空间 分 辩 率 的 遥感 影像 信息 不 再 只 由 一 
种 尺度 来 表示 ,可 由 多 种 适宜 的 尺度 来 描述 。 
在 多 尺度 分 割 中 ,由 于 尺度 因子 颜色 因子 、 形 状 因 
子 以 及 影像 各 波段 反射 率 权重 值 均 会 对 影像 分 割 的 
结果 产生 影响 ,尤其 是 尺度 因子 。 因 此 ,运用 最 
大 面积 法 “和 精度 对 比 法 ”确定 各 类 荒漠 地 物 的 
矿 度 因子 (最 优 分 割 太 度 ) 。 由 于 获取 的 影像 为 可 
见 光 ,在 3 个 波段 范围 内 ,地 物 都 有 相应 的 反射 特 
征 。 因 此 ,权重 值 均 为 1/3 。 而 颜色 因子 和 形状 因 
-F ,根据 实际 地 表 和 覆盖 物 颜 色 及 形状 , 设 定 为 0.6 和 


1.3.3 类 层次 结构 构建 ”根据 研究 区 地 表 分 类 体 
系 , 构 建 类 层次 。 在 轻 度 沙 漠 区 ,第 1 层次 ,分 市 尺 
度 为 15 ,利用 植被 指数 归 一 化 绿 红 差 异 指数 ( NGR- 
DI) 和 归 一 化 绿 蓝 差异 指数 (NGBDI) 区 分 植被 与 非 
植被 ;第 2 层次 ,对 第 1 层次 中 的 植被 进行 再 分 类 ， 
运用 植被 指数 .光谱 纹理 3 类 特征 ,在 植被 基础 上 提 
取 灌 木 ,未 分 类 的 为 草本 ;第 3 层次 ,对 第 1 层次 中 的 
非 植 被 进行 再 分 类 ,分 割 尺度 为 60 ,运用 光谱 和 纹理 
特征 提取 结 皮 ,未 分 类 为 裸 沙 。 中 度 沙 漠 化 区 和 重度 
沙漠 化 区 类 层次 结构 与 轻 度 沙漠 区 相似 (图 3)。 


表 1 影像 分 类 解 译 标志 


Tab.1 Interpretation symbols of image classification 


分 类 遥感 影像 特征 影像 显示 
"A 形状 规则 , 片 状 分 布 ,边界 清晰 , 呈 
土 黄色 。 部 分 地 区 出 现 风蚀 沟 餐 


形状 规则 ,条 带 状 分 布 ,边界 清晰 ， 


灌木 。 形状 较 规则 ,边界 清晰 , 插 深 绿 
O^ ERR, IH 


(S 


形状 不 规则 ,边界 不 清晰 , 呈 浅 绿 
[US D NASCE 


0.5, 
1.3.2 分 类 体系 ”在 前 人 研究 基础 上 ,通过 野外 调 
无 人 机 影像 
轻 度 沙漠 化 区 中 度 沙漠 化 区 重度 沙漠 化 区 
分 割 尺度 15, 分 割 尺度 15， 
植被 指数 (NGRDI NGBDI) GÉ 植被 指数 (NGRDI NGBDI) 


us 


图 3 多 尺度 分 割 分 层 


分 割 尺 度 20， 
植被 指数 (NGRDI. NGBDD 


分 割 尺度 70, 光 谱 植被 指数 
光 


非 植被 


分 割 尺 度 65, 光 谱 


特征 . 纹理 特征 谱 特 征 特征 .纹理 特征 


-:3IEZ 


Ds 


Fig.3  Layered chart of multi-scale segmentation 
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1.3.4 影像 对 象 特征 选取 及 信息 提取 采用 规则 
分 类 方法 对 影像 进行 分 类 ] 。 在 前 人 研究 基础 
EO ,根据 研究 区 情况 ,选取 以 下 10 个 指标 ( 表 2) 。 


R2 影像 对 象 特征 参数 


Tab.2 Characteristic parameters of image objects 


特征 类 型 影像 对 象 特征 特征 数量 /个 
光谱 特征 Mean( Red) 1 
Mean( Green) 1 
Mean( Blue) 1 
Brightnes 1 
Max. diff 1 
形状 特征 Length/Width 1 
Shape index 1 
纹理 特征 GLCM Homogeneity 1 
植被 指数 NGRDI 1 
NGBDI 1 

合计 10 


注 :Mean(Red) „Mean ( Green) , Mean( Blue) 分 别 为 对 象 在 波段 ( 红 、 
绿 、 蓝 ) 光 谱 的 均值 ;Brightness 为 对 象 亮度 值 ; Max. diff 为 最 大 化 差 
SE ERE; Length/ Width 为 对 象 的 长 宽 比 ;Shape index 为 形状 指标 ; GL- 
CM Homogeneity 24 JK BE3t/E Ali [Ve — 逆差 矩 ;NGRDI 为 归 一 化 绿 红 差 
异 指 数 ;NGBDI 为 归 一 化 绿 蓝 差异 指数 。 下 同 。 


1.3.5 ”结果 对 比 与 精度 评价 ”通过 混淆 和 矩阵 方法 
对 面向 对 象 多 尺度 分 割 信 息 提取 方法 与 基于 像 元 的 
最 大 似 然 法 信息 提取 方法 进行 结果 对 比分 析 5) 。 


2 结果 与 分 析 


2.1 地 表 类 型 的 分 割 尺 度 分 析 与 选取 

选择 10 ~100 之 间 的 尺度 对 影像 进行 分 割 , 以 
10 为 1 个 间隔 进行 9 次 分 割 实 验 (每 1 次 分 割 基于 
原始 数 进行 分 割 ) (图 4)。 在 9 次 分 割 实验 中 ,3 UE 


x p 
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段 反射 率 权 重 均 为 173 ,颜色 因子 和 形状 因子 分 别 
选择 0.6 和 0.5。 

由 图 4 可 见 ,分 制 尺 度 10 和 20 分割 对 象 较 破 
碎 , 斑 块 较 小 。 分 割 尺度 三 30 后 , 斑 块 较 大 , 均 出 
现 “ 过 分 割 ” 的 现象 。 表 明 分 割 尺 度 越 小 ,分 割 对 
象 斑 块 越 小 , 越 破碎 , 斑 块 数量 增多 ;分 割 尺 度 越 
大 ,分 制 对 象 斑 块 越 大 ,破碎 程度 小 , 斑 块 数量 
2.1.1 最 大 面积 法 “对 于 某 种 地 物 ,影像 分 割 后 该 
类 地 物 的 影像 对 象 都 存在 最 大 面积 .最 小 面积 及 平 
均 面 积 。 前 人 研究 ”1 表明 ,最 大 面积 与 分 割 尺 度 
有 一 定 的 规律 ,影像 分 割 的 对 象 不 应 大 于 实际 地 物 
的 大 小 ,能 充分 显示 实体 地 物 的 面积 、 空 间 结构 等 特 
点 。 最 优 尺度 不 是 一 个 点 ,而 是 一 个 区 间 范 围 。 

分 割 尺 度 设 定 在 10 ~100 间 , 以 5 为 1 个 区 间 ， 
统计 每 个 分 割 尺 度 下 影像 对 象 的 最 大 面积 。 

通过 图 5 可 知 , 随 着 分 割 尺度 的 增加 ,分割 斑 块 
最 大 面积 总 体 呈 上 升 趋势 。 在 图 Sa 20 ~ 50、 
60 ~75 .80 ~ 90 尺度 区 间 斑 块 最 大 面积 呈 上 升 态 
势 ; 在 15 ~20.55 -60,90 ~ 100 尺度 区 间 斑 块 最 大 
面积 较为 稳定 。 在 图 5b 中 20 ~70,80 ~90,95 ~ 
100 尺度 区 间 斑 块 最 大 面积 呈 上 升 态势 ;在 15 ~20、 
70 ~80 .90 ~95 尺度 区 间 斑 块 最 大 面积 较为 稳定 。 
在 图 5c 中 20 265,70 ~75 .80 ~85 区 间 尺 度 区 间 斑 
块 最 大 面积 呈 s 上 升 态势 ;在 13 ~20 .65 ~70、75 ~ 
80 .90 ~ 100 尺度 区 间 斑 块 最 大 面积 较为 稳定 。 表 
明 在 尺度 增加 时 ,分 割 对 象 在 减少 ,面积 呈 增 加 趋 
势 。 但 在 图 $a 中 75 ~80 区 间 和 图 5c 中 85 ~90 区 
H E FEBR, 由 于 多 尺度 分 割 是 根据 图 中 地 物 的 


图 4 不 同 尺 度 分 割 结果 图 


Fig.4 Results of segmentation on different scales 
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形状 颜色、 尺度 等 确定 不 同 地 物 的 分 割 边界 ,尺度 
增 大 时 分 割 对 象 的 数量 在 减少 ,但 分 制 对 象 的 最 大 
面积 不 一 定 会 增 大 。 


在 尺度 为 60 的 时 候 达 到 了 84. 32% ;在 图 6b 中 , 灌 
木 和 草本 在 尺度 20 分 割 的 精度 最 高 , 分别 为 
82.91% 和 82. 4396 , 结 皮 在 尺度 为 70 的 时 候 达 到 


根据 前 人 的 研究 汪汪 ,分 割 斑 块 面积 大 小 不 
超过 地 物 实体 大 小 , 即 在 一 定 尺 度 范 围 内 ,最 大 面积 
保持 不 变 , 该 类 地 物 的 实际 大 小 能 较 好 地 表现 。 选 
取 最 佳 分 割 尺 度 区 间 轻 度 沙 漠 区 在 15 ~ 20,55 ~ 
60,90 ~ 100; 中 度 沙漠 化 区 在 15 ~ 20,70 ~ 80, 
90 ~95 ;重度 沙漠 化 区 在 15 ~ 20,65 ~ 70,75 ~ 80, 
90 ~ 100, 

由 于 灌木 和 草本 个 体 较 小 ,尺度 较 大 易于 造成 
“从 分 割 " 现象 ; 裸 沙 和 结 皮 个 体 较 大 ,尺度 较 小 造 
成 “过 分 割 ” 的 现象 。 因 此 , 轻 度 沙漠 区 中 灌木 和 草 
本 分 割 尺 度 区 间 15 ~ 20, 结 皮 分 制 尺度 区 间 为 
55 ~ 60; 在 中 度 沙漠 化 区 中 灌木 和 草本 分 割 尺 度 区 
间 15 ~ 20 , 结 皮 分 割 尺度 区 间 为 70 ~ 80; 在 重度 沙 
漠 化 区 灌木 和 草本 分 割 尺 度 区 间 15 ~ 20 , 25 Rz 4 8l 
尺度 区 间 为 65 ~70 .75 ~80。 

2.1.2 精度 对 比 法 ”每 类 地 物 各 选择 80 个 地 物 样 
本 ,根据 尺度 区 间 ,对 比分 制 精度 。 

由 图 6 可 知 ,在 图 6a 中 ,灌木 和 草本 在 尺度 15 

分 割 的 精度 最 高 ,分 别 为 85.76% 和 83. 0496 , 结 


f 83.8196 ;在 图 6c 中 ,灌木 和 草本 在 尺度 15 分 制 
的 精度 最 高 ,分 别 为 84.05% 和 83. 30% , 结 皮 在 尺 
度 为 65 的 时 候 达到 了 84.87% 。 

因此 ,确定 轻 度 沙 漠 区 灌木 和 草本 分 割 尺 度 为 
15 , 结 皮 为 60; 中 度 沙漠 化 区 灌木 和 草本 分 割 尺 度 
为 20, 结 皮 为 70; 重 度 沙 漠 化 区 灌木 和 草本 分 割 尺 
度 为 15 , 结 皮 为 65 。 
2.2. ”地表 类 型 的 影像 特征 与 分 类 规则 
裸 沙 . 结 皮 灌木 .草本 4 类 地 物 分 别 选 取 100 
个 样本 ,统计 样本 的 光谱 形状 纹理 .植被 指数 4 类 
特征 值 ,对 特征 值 进行 统计 ,依据 95% 的 置信 区 间 
确定 阔 值 ,建立 规则 。 

从 表 3、 表 4 和 表 5 看 出 ,4 种 地 表 和 覆盖 物 的 光 
谱 和 形状 纹理 特征 相似 ,难以 区 分 。 植 被 指数 可 以 
准确 区 分 植被 与 非 植 被 ,但 灌木 与 草本 的 植被 指数 
相似 ,无 法 使 用 植被 指数 区 分 ,中 度 沙漠 化 区 和 重度 
沙漠 化 区 存在 相同 的 问题 。 单 一 特征 无 法 进行 准确 
区 分 ,因此 ,同时 选用 多 个 特征 进行 提取 ( 表 6、 表 7 
MKR 8)。 
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Fig.6 The accuracy of cl 
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图 6 不 同 分 割 尺度 的 分 类 精度 
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于 旱 区 研究 
A3 轻 度 沙漠 化 区 地 类 样本 特征 值 统计 
Tab.3 Statistical results of the eigenvalues of sample plots in slightly desertified area 

特征 值 梨 沙 (最 小 值 - 最 大 值 ) 结 皮 ( 最 小 值 -最 大 值 ) 灌木 (最 小 值 - 最 大 值 ) 草本 (最 小 值 - 最 大 值 ) 
Mean( Red) (71,131) (64,122) (26,99) (51,100) 
Mean( Green) (98,158) (78,137) (35,111) (61,117) 
Mean( Blue) (119,181) (88,151) (41,124) (69,134) 

Brightnes (96,154) 71,130 38,110 (61,115) 

Max. diff (0.2,0.5) (0.14,0.42) (0.14 ,0.6) (0.15,0.58) 
Length/ Width (1,4.26) (1.02,3.9) (1,4.32) (1,3.9) 
Shape index (1.01,2.45) (1.05,2.29) (1,2.71) (1.03,2.54) 

GLCM Homogeneity (0.009,0.5) (0.01,0.39) (0.005,0.17) (0.016,0.24) 
NGRDI ( -0.12, -0.04) ( -0.13, -0.02) ( -0.18, -0. 12) ( 20.14, -0.11) 
NGBDI (0.07 ,0.17) (0.03 ,0. 15) (0.19,0.25) (0.18,0.21) 

表 4 中 度 度 沙漠 化 区 地 类 样本 特征 值 统计 
Tab.4 Statistical results of the eigenvalues of sample plots in moderately desertified area 

特征 值 裸 沙 (最 小 值 - 最 大 值 ) 结 皮 ( 最 小 值 -最 大 值 ) 灌木 (最 小 值 - 最 大 值 ) 草本 (最 小 值 - 最 大 值 ) 
Mean( Red) (76 ,133 ) (69,125) (73,109) (60,144) 
Mean( Green) (94,159) (72,141) (76,129) (65,142) 
Mean( Blue) (109,188) (80,158) (85,155) (75,175) 

Brightness (95,158) (74,141) (80,128) (71,140) 

Max. diff (0. 14,0. 47) (0.12,0.46) (0.15,0.49) (0.17,0.56) 
Length/ Width (1,5.9) (1,4.8) (1.01,2.28) (1.01,3.96) 
Shape index (1.01,2.54) (1.01,2.35) (1.1,2.1) (1.0,2.29) 

GLCM Homogeneity (0.04,0. 08) (0.03 ,0.68) (0.04 ,0.58) (0.02 ,0.55) 
NGRDI (-0.03, -0.01) ( -0.03, -0.02) ( -0.13, -0.04) ( -0.11, -0.03) 
NGBDI (0.02,0.17) (0.02 ,0.2) (0.007 ,0.013) (0.01,0.02) 

表 5 重度 沙漠 化 区 地 类 样本 特征 值 统计 
Tab.5 Statistical results of the eigenvalues of sample plots in seriusly desertified area 

特征 值 梨 沙 (最 小 值 - 最 大 值 ) 结 皮 ( 最 小 值 -最 大 值 ) 灌木 (最 小 值 - 最 大 值 ) 草本 (最 小 值 - 最 大 值 ) 
Mean( Red) (64,124) (41,115) (30,92) (46,125) 
Mean( Green) (85,158) (56,143) (38,117) (62,156) 
Mean( Blue) (103,190) (68,171) (49,143) (83,187) 

Brightness (85,156) (56,143) (41,117) (64,155) 

Max. diff (0.31,0.53) (0.21,0.60) (0.22,0.75) (0.31,0.68) 
Length/ Width (1.00,8) (1.00,5.3) (1.00,5.65) (1.00,3.8) 
Shape index (1.00,5.2) (1,2.8) (1,2.8) (1.00,2.6) 

GLCM Homogeneity (0.17,0. 73) (0.13,0.65) (0.11,0.61 (0.16,0.66) 

NGRDI ( 20.13, -0.05) ( 20.14, -0.04) ( 20.18, -0.14) ( 20.16, -0.14) 

NGBDI (0.07,0. 16) (0.05 ,0.16) (0.24 ,0.37) (0.16,0.26) 

36 轻 度 沙漠 化 区 分 类 体系 与 规则 
Tab.6 Classification system and rules of slightly desertified area 
类 层次 ”分割 尺度 地 表 类 型 规则 
Level 1 15 植被 -0.18 <NGRDI< -0.12,0. 18 < NGBDI «0.25 
Level 2 15 灌木 26 < Mean( Red) < 95,38 < Brightness < 100, — 0. 18 < NGRDI < - 0. 12,0. 005 < GLCM Homogeneity < 0. 17, 
0.18 < NGBDI «0. 25 
Level 3 60 结 皮 78 «Mean(Green) «137,88 < Mean( Blue) «151 ,0. 14 < Max. diff «0.46,0.003 < GLCM Homogeneity «0.42 


201909.00074v1 


chinaXiv 


芯 佳 忆 等 :基于 无 人 


ChinaXiv 合 作 期 刊 


影像 的 荒漠 地 表 类 型 信息 提取 


R7 中 度 沙漠 化 区 分 类 体系 与 规则 


Tab.7 Classification system and rules of moderately desertified area 


类 层次 DARE 地 表 类 型 


规则 


Level 1 20 植被 -0.13 «NGRDI < -0.03,0.007 < NGBDI «0.22 
Level 2 20 灌木 ”76 < Mean( Green) «129,80 < Brightness < 128 ,0. 15 < Max. diff <0. 49 \0.007 < NGBDI «0.013 
Level 3 70 结 皮 69 < Mean(Red) «125,80 < Mean( Blue) «158,74 < Brightness < 141 ,0. 04 < GLCM Homogeneity « 0. 08 


表 8 重度 沙漠 化 区 分 类 体系 与 规则 


Tab.8 Classification system and rules of seriously desertified area 


类 层次 “分割 太 度 地 表 类 型 


规则 


Level 1 15 植被 


-0.18 <NCRDI< -0.13 ,0.16 <NCBDI<0.37 


Level 2 15 灌木 ”30 < Mean(Red) «92,38 < Mean(Green) «117,41 < Brightness < 117 ,0. 24 < NGBDI «0.37 


Level 63 SIE 0,13 <GLCM Homogeneity < 0. 65 

2.3 ”地 表 类 型 的 解 译 精 度 评价 与 对 比 
采用 以 上 分 类 方案 及 其 规则 ,分 类 结果 见 图 

7d ~7f, 图 8d ~8f 为 局 部 图 。 运 用 基于 像 元 的 最 大 
似 然 法 提取 地 表 类 型 信息 ,分 类 结果 见 图 7a ~7c， 
图 8a ~8c 为 局 部 图 。a 和 c 为 轻 度 沙漠 化 区 (固定 
Vb Fr.) ,b I c 为 中 度 沙 漠 化 区 ( 半 固 定 沙丘 ) ,d RIT 
为 重度 沙漠 化 区 (流动 沙丘 ) 。 每 类 地 物 均 随机 选 
择 150 个 地 物 ,采用 混淆 矩阵 ( 表 9、 表 10) 对 结 
行 精度 评价 与 对 比 。 


A EN ULM LAM v ws 
到 7 面向 对 象 的 多 尺度 分 割 方法 (d ~f) 和 基于 像 元 的 最 大 似 然 法 (a ~c) 信 息 提 取 结 果 


41 <Mean(Red) «115,68 < Mean( Blue) «171,56 < Brightness < 143 ,1. 00 < Length/Width «5.3, 


可 以 看 出 轻 度 沙 漠 区 总 体 分 类 精度 为 
93.00% , 较 基 于 像 元 的 最 大 似 然 法 提高 10. 34% ， 
Kappa 系数 达到 了 0. 89 , 较 基 于 像 元 的 最 大 似 然 法 
提高 0.09 ;中 度 沙 漠 化 区 总 体 分 类 精度 为 91. 8396 ， 
较 基 于 像 元 的 最 大 似 然 法 提高 13. 67% ,Kappa 系数 
达到 了 0. 84 , 较 基 于 像 元 的 最 大 似 然 法 提高 0. 12; 
重度 沙漠 化 区 总 体 分 类 精度 为 93. 50% , 较 基 于 像 
元 的 最 大 似 然 法 提高 12. 50% , Kappa 系数 达到 了 
0.89, 较 基 于 像 元 的 最 大 似 然 法 提高 0. 1 。 结 果 显 
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Fig.7 Object-oriented multi-scale segmentation and rule classification method (d-f) and maximum likelihood method based on 


pixel information extraction results (a-c) 
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图 8 面向 对 象 的 多 尺度 分 割 方法 (d ~f) 和 基于 像 元 的 最 大 似 然 法 (a ~c) 信 息 提取 结果 (局 部 图 ) 


Fig.8 Object-oriented multi-scale segmentation and rule classification method (d-f) and maximum likelihood method based on 


pixel information extraction results (a-c) (local diagram) 


干旱 区 研究 
RI 面向 对 象 的 多 尺度 分 割 方法 分 类 精度 评价 
Tab.9 Classification accuracy evaluation of object-oriented multi-scale segmentation method 
精度 类 别 轻 度 沙 漠 区 分 类 精度 中 度 沙 漠 化 区 分 类 精度 重度 沙漠 化 区 分 类 精度 

灌木 结 皮 裸 沙 草本 灌木 结 皮 裸 沙 草本 灌木 结 皮 裸 沙 草本 

灌木 140 1 3 7 138 4 0 9 142 1 3 J: 

结 皮 2 138 4 0 4 136 8 5 2 141 4 0 

裸 沙 0 8 143 6 0 4 142 0 0 5 143 4 

草本 8 3 0 137 8 6 0 136 6 3 0 139 

生产 者 精度 /% 93.33 92 95.33 91.33 92 90.67 94.67 90 92 95.33 95.33 91.33 
用 户 精度 /% 86 91.33 96.67 82 91.39 86.62 97.26 86.62 90. 79 91.08 99.3] 87.26 


总 体 精 度 =93. 00% 
Kappa 系数 =0. 89 


表 10 


总 体 精度 =91. 83% 


总 体 精度 =93.50% 


Kappa 系数 =0.84 Kappa 系数 =0. 89 


基于 像 元 的 最 大 似 然 法 分 类 精度 评价 


Tab.10 Classification accuracy evaluation of maximum likelihood method 


— 轻 度 沙 漠 区 分 类 精度 中 度 沙漠 化 区 分 类 精度 重度 沙漠 化 区 分 类 精度 

灌木 。 结 皮 AWV 草本 灌木 结 皮 RWV MAK 灌木 å 结 皮 AO 草本 

灌木 125 4 0 22 116 11 0 22 112 12 0 10 

结 皮 1 118 14 6 7 113 16 17 13 123 8 26 

裸 沙 3 11 136 5 6 15 130 1 4 2 142 5 

草本 21 17 0 117 21 11 4 110 21 13 0 109 
生产 者 精度 /% —— 82.78 84.89 87.74 75.48 77.85 73.86 85.53 75.34 83.58 72.35 92.81 76.22 
用 户 精度 /% 83.33 78.67 90.67 78 77.33 75.33 86.67 73.33 74.67 82 94.6] 72.61 


总 体 精 度 — 82. 66% 
Kappa 系数 =0. 80 


示 本 方法 分 类 精度 高 于 基于 像 元 的 最 大 似 然 方法 ， 
表明 本 方法 能 较为 精确 地 提取 荒漠 地 表 类 型 信息 ， 
在 高 分 辩 率 遥感 影像 荒漠化 信息 提取 中 具有 一 定 的 
优势 。 


3 讨论 与 结论 


3.1 讨论 

该 方法 分 类 效果 受到 “影像 分 制 "” 和 “特征 提 
取 ” 两 个 过 程 影响 较 大 。 分 割 尺度 决定 影像 对 象 的 
大 小 及 信息 提取 的 精度 ,而 尺度 没有 明确 的 单位 , 值 
越 大 则 分 割 的 影像 对 象 的 斑 块 越 大 ,反之 则 越 小 。 
采用 最 大 面积 法 和 精度 对 比 法 5 2 确定 分 割 尺 度 ， 
在 一 定 程度 上 能 避免 “过 分 制 " 与 “ 欠 分 制 " 现象。 
由 于 高 分 遥感 影像 具有 “ 同 谱 异 物 , 同 物 异 谱 ” 的 特 
性 ,因此 ,规则 的 建立 能 够 有 效 地 区 分 地 物 , 提 高 分 
类 精度 ,但 不 同 荒漠 化 程度 中 地 表 类 型 特征 不 同 , 需 
选用 多 种 特征 组 合 对 地 物 进 行 提取 。 在 前 人 研究 基 
MEO? ,选取 光谱 ,形状 ,纹理 ,植被 指数 4 种 特 
征 建立 规则 ,有 效 地 避免 此 现象 出 现 ,分 类 精度 得 到 
了 明显 提高 ,这 与 黄 慧 萍 2 的 研究 一 致 。 但 错 分 现 
象 还 是 存在 ,如 能 结合 地 物 分布 区 域 的 地 形 数据 , 完 
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总 体 精 度 =78. 16% 


总 体 精度 =81.00% 


Kappa 系数 =0.72 kappa 系数 =0.79 


善 典型 地 物 影 像 特 征 库 ,选取 具有 客观 性 和 代表 性 
的 特征 ,分 类 效果 应 该 会 更 好 ,精度 会 更 高 。 

面向 对 象 的 多 尺度 分 割 方 法 分 类 精度 高 的 原 
:J 运用 高 分 辩 率 影像 特点 , 极 大 改善 了 异物 同 
谱 的 问题 ,有 效 避 免 了 椒盐 效应 ;@) 利用 光谱 纹理 
及 形状 特征 较为 准确 地 对 研究 区 不 同 地 物 进 行 描 
述 ,使 易 混淆 区 域 的 识别 度 更 高 。 综 上 所 述 ,在 最 优 
分 割 情况 下 地 物 特征 能 够 得 到 准确 描述 ,不 同 地 物 
之 间 的 差异 性 得 到 充分 体现 ,降低 影像 中 特征 相似 
地 物 之 间 的 误 分 现象 。 选 取 合 理 的 分 类 特征 ,建立 
分 类 规则 ,能 够 实现 信息 高 精度 提取 。 本 文 在 分 类 
体系 中 存在 一 定 的 瑕 羔 , 先 统 的 将 地 表 类 型 划 为 4 
类 ,不 能 准确 地 说 明芳 漠 化 演变 的 过 程 。 如 能 将 地 
表 类 型 进一步 细 分 ,分 析 地 表 类 型 如 何 变 动 , 可 进 一 


步 提高 精度 。 
3.2 结论 


基于 无 人 机 影像 对 古 尔 班 通 古 特 沙漠 地 表 类 型 
信息 进行 提取 ,运用 面向 对 象 的 多 尺度 分 割 ,结合 最 
大 面积 法 和 精度 对 比 法 进行 分 割 实验 ,确定 最 优 分 
割 太 度 。 在 此 基础 上 ,综合 样本 提取 的 光谱 形状、 
纹理 ,植被 指数 特征 建立 规则 ,并 采用 规则 分 类 的 信 
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息 提 取 地 表 类 型 信息 。 


(1) 荡 江 地表 类 型 不 同 ,最 佳 分 割 尺度 不 同 。 


小 尺度 适用 于 灌木 .草本 小 斑 块 荒漠 地 表 类 型 ;大 尺 


度 适 


用 于 裸 沙 i EK BEBO ADIT 
(2) 不 同 亮 漠 化 程度 地 表 类 型 特征 相似 ,植被 


指数 能 准确 区 分 植被 与 非 植被 ,而 裸 沙 、 结 皮 灌木 、 
草本 4 类 地 物 的 光谱 和 形状 、 纹 理 特 征 相 似 ,无 法 运 
用 单个 特征 进行 区 分 , 需 选用 多 种 特征 组 合 提取 地 
表 类 型 。 


(3) 本 方法 分 类 精度 高 于 基于 像 元 的 最 大 似 然 


法 ,表明 针对 可 见 光 无 人 机 影像 ,面向 对 象 多 尺度 分 


LI 


割 方法 能 高 精度 地 提取 荒漠 地 表 类 型 信息 。 
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Information Extraction of Desert Surface Types Based on UAV Image 
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Abstract: Based on the images of unmanned aerial vehicle ( UAV ) ,in this study the surface type information of 
the Gurbantunggut Desert was extracted , and the object-oriented multi-scale segmentation was used to extract the in- 
formation of the sample plots of surface types from the spectrum , shape, texture and vegetation index. The results 
showed that :CD The best segmentation scale for the different desert surface types was different; (2) The features of 
surface types with different desertification levels were similar and could not be distinguished by single characteris- 
tics;(3) Compared with the pixel-based maximum likelihood method, the object-oriented multi-scale segmentation 
method was improved significantly. The overall classification accuracies of the object-oriented slightly , moderately 
and seriously desertified areas were 93. 0096 ,91. 8396 and 93. 5096 respectively , and they were 10. 3496 ,11. 8696 
and 12.5046 higher than those of the pixel-based maximum likelihood method. The results revealed that the object- 
oriented multi-scale segmentation method could be used to extract the desert surface type information with high ac- 
curacy for the UAV visible image. 


Key words; UAV image; desert; land surface type; information extraction; the Gurbantonggut desert 


